Origin and roles of a strong electron-phonon interaction 
in cuprate oxide superconductors 
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A strong electron-phonon interaction arises from the modulation of the superexchange interac- 
tion by phonons. It is responsible for the softening of the half breathing modes, as is studied in 
Phys. Rev. B 70, 184514 (2004). When antiferromagnetic spin fluctuations are developed around Q, 
Cu-0 bond stretching modes can also be soft around twice Q due to the electron-phonon interaction. 
However, it can play no or only a minor role in the binding of d7-wave Cooper pairs. 

PACS numbers: 71.38.-k, 74.20.-z, 75.30.Et 
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It is an important issue to elucidate high critical 
temperature (high- Tc) superconductivity in cuprate ox- 
ides discovered in 1986 0. It occurs in the vicin- 
ity of the Mott metal-insulator transition or crossover. 
In 1987, a theory is proposed that Gutzwiller's quasi- 
particles are bound into d7-wave Cooper pairs due to 
the superexchange interaction between nearest neighbors 
0, According to observed specific heat coefficients 
7 ~ 14 mJ/K^mol 0, the bandwidth of quasiparticles 
is W* — 0.3-0.4 eV. Since the superexchange interaction 
constant is as strong as J = —(0.10-0.15) eV 0], ob- 
served high Tc can be easily reproduced. High- Tc super- 
conductivity occurs in an intermediate-coupling regime 
I J|/M^* = 0.3-0.5 for superconductivity, which is reahzed 
in the strong-coupling regime for electron correlations. 

It is also an important issue to elucidate the origin 
and roles of a strong electron-phonon interaction, whose 
presence is implied by pieces of evidence: the softening 
of the so called half breathing modes around (±7r/a, 0) 
and (0, ±7r/a) in the two-dimensional (2D) Brillouin zone 
(BZ) 0,01111113' with a the lattice constant of Cu02 
planes, the softening of Cu-0 bond stretching modes 
around (±7r/2a,0) and (Cl,±7r/2a) kinks in the 

quasiparticle dispersion |l3l [ij] , and so on. On the other 
hand, observed small isotope shifts of Tc imply that 
the strong electron-phonon interaction plays only a minor 
role in the occurrence of superconductivity. 

Since charge fluctuations are suppressed by strong elec- 
tron correlations, the conventional electron-phonon in- 
teraction arising from charge-channel interactions must 
be weak in cuprate oxides. An electron-phonon interac- 
tion arising from spin-channel interactions can be strong 
when electron correlations are strong. For example, it 
plays a significant role in the spin-Peierls effect. Then, a 
novel one is proposed in a previous paper [l6l | : one arising 
from the modulation of the superexchange interaction by 
phonons. It can explain the softening of breathing modes 
around (±7r/a,0) and (0, ±7r/a); it is predicted that the 
softening must be small around (±7r/a, ±7r/a). An at- 
tractive interaction arising from the virtual exchange of a 
single phonon of the modes is very weak between nearest 
neighbors; it is strong between next-nearest neighbors. 



Then, it can play no significant role in the formation of 
1^7- wave Cooper pairs |l6| . The observed small isotope 
shifts of Tc can never contradict the presence of this novel 
and strong electron-phonon interaction. One of the pur- 
poses of this Letter is to explain the softening around 
(±7r/2a,0) and (0,±7r/2a) in 2D-BZ. 

The superexchange interaction arises from the vir- 
tual exchange of pair excitations of 3d electrons across 
the upper Hubbard band (UHB) and the lower Hub- 
bard band (LHB) [13. When their non-zero bandwidths 
are ignored, the exchange interaction constant between 
nearest-neighbor Cu ions is given by 
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matrix between nearest 
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neighbor 2p and 3d orbits, and £p the depths of 3c? and 
2p levels, and U the on-site repulsion between 3c? elec- 
trons. One may argue that parent cuprate oxides with 
no hole doping must be charge-transfer insulators rather 
than Mott insulators because doped holes reside mainly 
at O ions, so that ep > e^. However, this argument dis- 
agrees with Erf — Ep ~ 1 eV predicted by band calculations 
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The preferential doping does not necessarily 
mean Cp > e^j but it simply means that the local charge 
susceptibility of 3c? electrons is much smaller than that of 
2p electrons. When we follow the band calculations and 
we use V=1.6 eV, ed-ep=l eV, and U=5 eV [l^, Eq. ((TJ 
gives J — —0.27 eV. This is about twice as large as the 
experimental J = —(0.10-0.15) eV 0. This discrepancy 
is resolved when nonzero bandwidths of UHB and LHB 
are considered [T^ . 

Displacements of the ith Cu ion and the [ij]th O ion, 
which lies between the ith and jth Cu ions, are given by 
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with Rj and R[y] = (l/2)(R.i-|-Rj) positions of the ith 
Cu and [zjjth O ions. Aid and Mp masses of Cu and O 
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ions, ^Aq and bx_q annihilation and creation operators 
of a phonon with a polarization A and a wave vector q 
or -q, ojxq a phonon energy, eAq = (cAq.x, ^Aq.^, CAq^z) 
a polarization vector, and N the number of unit cells. 
Here, we consider only longitudinal phonons; we assume 
EAq = (fe, 9y, 9z)/|q| for q within the first BZ. The q 
dependence of Wd.Aq and Wp.Aq is crucial. For example, 
'^d.Aq = and Wp,Aq = 0(1) for modes that bring no 
change in adjacent Cu-Cu distances. 

Denoting creation and annihilation operators of 3d 
electrons at the ith site by and die and those of 
3d electrons with wave number k by d^^^ and d^a, we 
define spin operators by Si = J2af3 h^'^^'^la'^'^P ^^'^ 

S(q) = (I/VN) Eko;3 5'^"''rf;k+iq)o^(k-iq)/3' ^ith CT = 

(ca; , CTj, , (Jz ) the Pauli matrixes. Two types of electron- 
phonon interactions arise from the modulations of J by 
the vibrations of O and Cu ions: 



Hp — iCp ^ ^ 
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X J2 ^r(q)7'r(q), 

r=s,d 



^Lq + ^Aq 
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with Cp and Cd given in the previous paper |l6|| . 

^s(q) = 2 [exq^x sin {qxa/2) + ex^^y sin (qya/2)] , (6) 

%(q) = 2 [exq^x sin {qxa/2) - exq^y sin {qya/2)] , (7) 

77s(k) = cos{kxa) + cos{kya), (8) 

r]d(k) = cos{kxa) — cos{kya), (9) 

being form factors, and 

^r(q) = ^ E '^Hq') [S (q'+iq) • S (-q'+iq)] , (10) 
q' 

being a dual-spin operator. Here, the x and y axes are 
within Cu02 planes, and the z axis is perpendicular to 
Cu02 planes. The d-p model is approximately mapped 
to the t-J model jlJ]. Then, we consider the i- J- infinite 
U model including Tip and Tid on a square lattice: 
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+ Uoo E 4^di^4l<^il + + Wd, (11) 
■i 

with (ij) over nearest neighbors, and Uoo an in- 
finitely large on-site repulsion to exclude double occu- 
pancy. Quasi-2D features are considered in terms of an 
anisotropy factor of spin fluctuations introduced below. 



Every physical quantity is divided into single-site and 
multi-site ones. Calculating the single-site one is reduced 
to determining and solving selfconsistently the Anderson 
model, which is an effective Hamiltonian for the Kondo 
problem. This is the single-site approximation (SSA) 
that includes all the single-site terms the 
SSA can also be formulated as the dynamical-mean-field 
theory 1251 or the dynamical-coherent-potential approx- 
imation [2g. Multi-site or intersite effects can be per- 
turbatively considered starting from a unperturbed state 
constructed in the non-perturbative SSA theory. Such a 
perturbative theory is nothing but a Kondo-lattice the- 
ory. 

The irreducible polarization function in spin channels 
is the sum of single-site TTs{iuji), which is the same as 
that of the Anderson model, and multi-site AtTj (iuji , q) : 
Trs{iuJi,q) = TTsiiuJi) -I- AtTs (iw; , q) . Spin susceptibilities 
of the Anderson model and the t-J model are given by 
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with J(q) = 2j77s(q). A physical picture for Kondo lat- 
tices is that local spin fluctuations at different sites inter- 
act with each other by an intersite exchange interaction. 
Following this picture, we deflne an intersite exchange 
interaction 1^ (iuii , q) by 
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Then, it follows that Is{iuJi,q) = J{q)-\-2U^ATTs{iuJi,q). 
The main term of 2U^ATrs{iuJi,q) is an exchange in- 
teraction arising from the exchange of pair excitations 
of Gutzwiller's quasiparticles jl^. It has a novel prop- 
erty that its strength is proportional to the width of 
Gutzwiller's band. It is antiferromagnetic (AF) when the 
chemical potential lies around the center of Gutzwiller's 
band or the nesting of the Fermi surface is sharp enough; 
it is ferromagnetic when the chemical potential lies 
around the top or bottom of Gutzwiller's band. 

In cuprate oxides, the exchange interaction /s(0, q) is 
AF and is maximal around nesting wave numbers of the 
Fermi surface. We assume that it is maximal at Q = 
(±37r/4a, ±n/a) and (±7r/a, ±37r/4a) in 2D-BZ, and that 
the susceptibility is approximately given by 



Xs (iw;,Q + q) = 



Xs(0,Q)«;2 



'■ + {q\\ay + 6'iq,cy + \u;i\/TAi 



(15) 



around each of Q, where q|| = {qx,qy), ^af is an energy 
scale of AF spin fluctuations, and c is the lattice constant 
along the z axis. The anisotropy factor 6 is introduced to 
consider quasi-2D AF spin fluctuations; the correlation 
length within the x-y plane is a/n and that along the 
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z axis is i5c//t. A cut-ofF qc = 7r/3a is introduced so 
that Xs ii^u Q + q) = for \qx\ > or \qy\ > q^. The 
anisotropy of the lattice constants plays no role when 5 
and qc are defined in these ways. 

The Green function for phonons is given by 
Dx{iuJi,q) = 2cjAq/[(jt^0^-^Aq + 2^AqAa;Aq(«t^i)], with 
AcjAq(iw;) the renormalization or softening of phonons. 
Three types of fluctuations can be developed: AF spin 
ones due to Is{iuji,q), and d'y-wave superconducting (SC) 
and charge bond order (CBO) ones due to a mutual in- 
teraction between quasiparticles, which is given by 

Xs[iuJi,q) - XsiiUJi) ^ - — TTT n~ (16) 

l-jIs{lUJi,Cl)Xs(lUJl) 

multiplied by and two single-site reducible vertex 
functions. The softening around 2Q occurs mainly due 
to AF spin ones around Q; 2Q's are equivalent to 
(±7r/2a, 0) and (0, ±7r/2a). When only the AF spin ones 
are considered, it follows that [l6| 
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In Eq. lUni), two Xs's appear because of the dual-spin 
operator (|10|) . Since we are interested in Cu-0 bond 
stretching modes around 2Q, we ignore vibrations of Cu 
ions; we assume that \CdVd,Xq}\/ Mp/Md — and 



\CpVp,xq\ = Cp eV/A, 



(20) 



where is a dimensionless constant and it is likely Cp — 
0(1) 16] . Since the contribution from small p is large 
in Eq. H19|) . we consider only the contribution from the 
r = s channel. The softening is given by 



(21) 



with 
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P \t*\ 



[Xs(0,Q)^'\t*\f meV, (22) 



E{iuJi,q) = ri^{q)v't (iq) 



^AF [Xs{Q)^'Y 



(23) 



In Eq. (|22|l . t* is the efi^ective transfer integral between 
nearest neighbors for quasiparticles, and we assume |i* | — 
W*/8 ~ 40-50 meV and WAq = 50 meV. It is likely that 
TAF/\t*\ = 0(1) and XsiO,Q)K^\t*\ = 0(1). We note 
that S(z<x';,q) is defined as a dimensionless quantity. 

We examine the k, i5, and to dependence of E{uj + iO, q); 
T = K is assumed in the oji' sum of Eq. p9|l . and 
the softening around one of 2Q's or 2Qo, with Qo = 
(— 37r/4a, tt/o) in 2D-BZ, is considered; 2Qo is equivalent 
to (7r/2a, 0). Figurenshows Re[S(w+i0, q)] as a function 
of or qx for several sets of S and lu/Taf- According 
to Fig.n Re[AtJAq(t^ + iO)] has a minimum around 2Qo 
as a function of q; it is likely that its minimum value is 
as large as -(10-20) meV for <^ 1 and (5 < 1. 

It is definite that ^ 1 in the critical region, and it 
is certain that S < 1Q~^ for cuprate oxides. Then, we 
propose that this second-harmonic effect of AF spin fiuc- 
tuations is mainly responsible for the observed softening 
[Ullii as large as -(10-20) meV around 2Q. 

Since the softening is small when is large or AF 
spin fluctuations are not developed, it must be small 
in the so called over-doped cuprate oxides. When AF 
spin fluctuations are developed similarly or differently 
between (±37r/4a, ±7r/a) and (±7r/a, ±37r/4a) because 
of the anisotropy of the Fermi surface within 2D-BZ, 
the softening must also occurs similarly or differently be- 
tween (±7r/2a,0) and (0,±7r/2a). These two predictions 
are consistent with observations [ill . . 

The soft modes with 2Q correspond to the so called 
stripes and checker-boards observed at low temperatures 
[UliiHillli- Since a charge density wave (CDW) 
can appear following the complete softening, it should 
be examined if the stripes and checker-boards can be ex- 
plained in terms of the CDW. In general, the 2Q com- 
ponent of the density of states, P2q(£), as a function of 
e is composed of symmetric and asymmetric ones with 
respect to the chemical potential or e = 0. The asym- 
metric one is large when CDW with 2Q is stabilized as 
a fundamental 2Q effect. According to an experiment 
|28|. the symmetric one is larger than the asymmetric 
one. This contradicts the scenario of CDW, even if the 
softening of the 2Q modes is large and the 2Q fluctu- 
ations are well developed; the amplitude of CDW must 
be small, even if it is stabilized. On the other hand, 
the symmetric one is large when the 2Q modulation is 
due to a second-harmonic effect of a spin-density wave 
(SDW) with Q (33. The observed almost symmetric 
P2Q (e) can be explained by the second-harmonic effect of 
SDW. When stripes and checker-boards are really static 
orders, stripes must be due to single-Q SDW and checker- 
boards must be due to double-Q SDW; magnetizations 
of the two waves must be orthogonal to each other in 
double-Q SDW 

On the other hand, it is proposed that a stripe or 
a checker-board at rather high temperatures must be an 
exotic ordered state, that is, a fluctuating state in a quan- 
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(a) (iv) (b) (iv) 




{q,+2nla)l{nla) 



FIG. 1; (a) Re[^(tj + jO,2Qo)] as a function of , with Qo = (— 37r/4a, 7r/a), and (b) Re[H(tj + iO, q)] as a function of 

(-27r/a < (/^ < -vr/a) for qy = 2-R/a. For the anisotropy factor, (i) 5=1, (ii) S = 10"^/^, (iii) S = 10"\ and (iv) 5 = 10"^. In 
each figure, solid, dotted, dashed, and dashed chain lines are for uj/Taf ~ 0.2, 0.4, 0.8, and 1.6, respectively. 



turn disordered phase. It should be examined whether it 
is actually such an exotic state. Another possibility is 
that it is a rather normal low-energy fluctuating state, 
whose energy scale is as small as that of the soft phonons. 
The other one is that it is an disorder- induced SDW ; 
it must be a rather simple but inhomogeneous SDW. 



According to Eq. (|I6|I . two mechanisms of attractive 
interactions, the spin- fluctuation one and the exchange- 
interaction one, are essentiaUy the same as each other in 
Kondo lattices. However, the main part of the attractive 
interaction in cuprate oxides is the superexchange inter- 
action rather than an interaction mediated by low-energy 
AF spin fluctuations. Since it is as strong as J = —(0.10- 
0.15) eV 0, observed high Tc can be easily reproduced. 
Although the strong electron-phonon interaction plays 
only a minor role in the formation of d'y-wave Cooper 
pairs, a small isotope shift of Tc can arise from the de- 
pression of superconductivity by the 2Q fluctuations. 



In conclusion, the strong electron-phonon interaction 
arising from the modulation of the superexchange inter- 
action by phonons is strong in cuprate oxides supercon- 
ductors. It is responsible for the softening of the half 
breathing modes around (±7r/a,0) and (0,±7r/a). In the 
critical region of SDW, where antiferromagnetic spin fluc- 
tuations are developed around nesting wave numbers Q 
of the Fermi surface, Cu-0 bond stretching modes can 
also be soft around 2Q. The softening is accompanied 
by the development of 2Q or 2Q x 2Q fluctuations, that 
is, stripe or checker-board fluctuations. However, the ob- 
servation that in ordered stripe and checker-board states 
the 2Q component of the density of states is almost sym- 
metric with respect to the chemical potential can never 
be explained by CDW with 2Q following the complete 
softening of the 2Q modes; they can be explained by the 
second- harmonic effect of SDW with Q. 
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